The Adsorption of Atomic Nitrogen on Ru(OOOl): Geometry and Energetics 
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The local adsorption geometries of the (2 x 2)-N and the (\/3 x \/3)I^30°-N phases on the 
Ru (0001) surface are determined by analyzing low-energy electron diffraction (LEED) intensity 
data. For both phases, nitrogen occupies the threefold hop site. The nitrogen sinks deeply into the 
top Ru layer resulting in a N-Ru interlayer distance of 1.05 A and 1.10 A in the (2 x 2) and the 
(\/3 X \/3)R30° unit cell, respectively. This result is attributed to a strong N binding to the Ru 
surface (Ru-N bond length = 1.93 A) in both phases as also evidenced by ab-initio calculations 
which revealed binding energies of 5.82 eV and 5.59 eV, respectively. 
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The rate-limiting step in the Haber-Bosch synthesis 
(ammonia synthesis) is represented by the dissociation 
process of molecular nitrogen over the active metal sur- 
faces [1]. While for the Fe(lll) surface the dissociative 
sticking coefficient of nitrogen is about 10"^ [2], the Ru 
surfaces exhibit only an exceedingly small sticking co- 
efficient of about 10^^^ at room temperature [3]. Yet, 
the promoted Ru catalyst is even more efficient than the 
industrially used Fe catalyst [4] since blocking of the sur- 
face by too much nitrogen and/or poisoning by traces of 
oxygen is less severe. Apart from the low sticking prob- 
ability, the turn over rates of the ammonia reaction are 
determined by the energy barriers of adsorbed hydrogen 
and nitrogen to form NH^;. This barrier is governed by 
the adsorption energy and the adsorption geometry of 
atomic nitrogen since N is much more strongly bound 
to Ru than H. In this Letter, we present the atomic ge- 
ometry and energetics of the atomic nitrogen adsorbed 
on Ru (0001) at two distinct coverages by employing the 
techniques of low-energy electron diffi'action (LEED) and 
density- functional theory (DPT). 

The experiments were conducted in a UHV system 
(base pressure 2.0 x 10"'^° mbar) equipped with a 4- 
grid LEED optics, a CMA for Auger spectroscopy and 
a quadrupole mass spectrometer for thermal desorp- 
tion spectroscopy (TDS). The Ru (0001) crystal could be 
heated resistively up to 1500 K and cooled to 130 K us- 
ing liquid nitrogen; the temperature was measured with 
a NiCr-Ni thermocouple spot welded to the back of the 
sample. Further details about the experimental setup 
can be found elsewhere [7]. The Ru(OOOl) crystal was 
cleaned in a standard way by Ar+ sputtering and re- 
peated cycles of oxygen exposure followed by annealing to 
1560 K in order to remove carbon contamination. Since 
the dissociative sticking coefficient of molecular nitrogen 
is extremely small we used ammonia and its ability to 
decompose at the Ru(OOOl) surface above 350 K into 
atomic nitrogen and hydrogen which desorbs at the cho- 
sen sample temperature. Nevertheless, NH3 exposures of 
several thousands of Langmuirs were needed to accom- 
modate substantial amounts of atomic nitrogen at the 



Ru surface. 

More specifically, the atomic nitrogen overlayers were 
prepared by adsorbing ammonia, which is known to disso- 
ciate into N and NH already during chemisorption above 
the desorption temperature of NH3. There is a delicate 
balance for the uptake of atomic nitrogen: NH3 dos- 
ing at lower temperatures results in a higher saturation 
coverage of atomic nitrogen but at the expense of very 
high NH3 doses, while at higher adsorption temperature 
the uptake of atomic nitrogen takes place rapidly reach- 
ing a smaller saturation coverage. To get the required 
high N coverage on a reasonable time scale, the temper- 
ature was lowered step by step starting from 580 K to 
room temperature while dosing a total amount of about 
20000 L of NH3 by backfilling the UHV chamber with 
1.0 x 10~^ mbar NH3. The specific protocol of NH3 
dosing can be found in Ref. [8]. After this procedure 
it took typically 12 hours to retain a base pressure of 
2.0 X 10"^ mbar (the main constituent of the residual 
gas was still due to NH3) as ammonia easily sticks to 
and desorbs from the walls of the UHV chamber. 

The (V3 X V3)R30° and (2 x 2) phases of nitro gen were 
prepared by dosing NH3 as described above and by heat- 
ing the sample to 525 K and 615 K for the (\/3 x V3)R30° 
and the (2 x 2) phase, respectively (heating rate 5 K/s). 
From HREELS measurements it is known that for tem- 
peratures beyond 500 K neither NH nor NH2 species are 
left on the surface [8]. 

After coohng the sample to 140 K LEED exhib- 
ited well-ordered (2 x 2) and (VS x y3)R30° patterns. 
LEED I(V) curves were taken at normal incidence. Pour 
fractional-order beams and two integer-order beams in an 
energy range from 50 eV to 300 eV were collected in the 
case of the (2 x 2), three fractional-order beams and two 
integer-order beams in an energy range from 30 eV to 
350 eV for the (VS x -\/3)R30°. The measurements were 
recorded on-line from the fluorescent screen of a display 
type 4-grid LEED optics by using a computer-controlled 
video-LEED system [9] . The experimental data were fur- 
ther processed by averaging over symmetrically equiva- 
lent beams, smoothing, and correcting for constant emis- 
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sion current. 

Due to re-adsorption of NH3 from the residual gas 
it was necessary to heat the sample again to 525 K 
[( V3 X V3)R30°-N] or 615 K [(2 X 2)-N] after a cooling and 
measuring period of about 5 minutes. This precaution 
ensures that the NHs-induced changes in the LEED I(V) 
curves of fractional-order beams remained small as quan- 
tified by a Pendry r-factor [10] of Rp < 0.10; note that 
the experimental uncertainty due to re-preparation was 
found to be about Rp w 0.08. 

To evaluate the experimental I(V) data, fully- 
dynamical LEED calculations were performed using the 
program package of Moritz [11]. Nine spin-averaged, rel- 
ativistically corrected phase shifts were used for the Ru 
substrate [12], and the phase shifts of nitrogen were taken 
from the literature [13] which have been already success- 
fully applied to the analysis of molecular nitrogen on the 
clean and 0-precovered Ru (0001) surface [14]. The influ- 
ence of thermal vibrations on the LEED intensities was 
accounted for by correcting the phase shifts with a fixed 
effective Debye temperature of 400 K for the substrate 
and for nitrogen. 

The agreement between experimental and calculated 
I(V) curves was quantified by the RDE factor introduced 
by Kleinle et al. [15] and by Pendry's r-factor Rp [10]. 
The refinement of structural parameters and of the real 
part of the muflan-tin zero was carried out by employing 
a nonlinear least-squares optimization scheme [16] based 
on the least-squares sum either of the LEED intensities 
or of the corresponding Y functions. 

For the density-functional theory (DFT) calculations 
we employed the generalized gradient approximation 
(GGA) of Perdew et al. [17] for the exchange-correlation 
functional and used ab-initio pseudo-potentials created 
by the scheme of TrouUier and Martins in the fully sep- 
arable form [18]; the electronic wave functions were ex- 
panded in a plane- wave basis set. The surface was mod- 
eled using the supercell approach, incorporating five lay- 
ers of Ru (0001) where the N atoms were adsorbed on 
one side of this slab [19]. Calculations were performed 
using a (2 X 2), a (V3 X V3)R30° and a (1x1) surface 
unit cell with an energy cut-off of 50 Ry. The integral 
over the Brillouin zone was performed using a special k- 
point set [20] , and care was taken to ensure an equivalent 
sampling in all (surface) geometries, with 9 k-points in 
the irreducible part of the (1x1) Brillouin zone. Fermi 
broadening of the occupation numbers was done with a 
width of 0.1 eV. The calculation scheme allows for simul- 
taneous relaxation of the electrons and atoms, where we 
relaxed the positions of the N atoms and the atoms in the 
top Ru layer, keeping the lower three Ru layer spacings 
fixed at the bulk values. 

To estimate the accuracy of the calculations, we 
studied the clean substrate Ru (0001) and the free N2 
molecule. The lattice constant of Ru obtained turned out 
to be 2.78 A, slightly larger than in experiments (2.71 A) 



as usually observed with GGA. Therefore the c/a ratio 
was somewhat smaller than the experimental value, be- 
ing 1.5755. When relaxing the (0001) slab, the outermost 
layer spacing contracts by about 3 %, a value which is 
also slightly larger than the experimental value of 2±1 % 
[21]. For the very strong bond of the N2 molecule we 
obtained a binding energy of 10.24 eV (not corrected for 
zero-point motion), a N-N bond length of 1.104 A, and 
a N-N vibrational frequency of 2336 cm~^. These values 
compare well with the experiment of values of: 9.77 eV, 
1.097 A, and 2359 cm"!. 

The LEED analyses were performed by placing the 
nitrogen in both phases in high-symmetry sites - i.e., 
bridge, on-top, threefold hollow fee and hep - and vary- 
ing the N-Ru layer spacing, the topmost Ru-Ru layer 
distance as well as the lateral and perpendicular displace- 
ments in the top Ru layer as long as they are compatible 
with the symmetry of the unit cell. To explore a wide 
range of parameter space, we changed, in a first step, 
the principal structure parameter, i. e., the N-Ru layer 
spacing, from 0.5 A to 2.5 A in steps of 0.1 A, keeping 
the other co-ordinates of the Ru atoms at those values 
found for the clean Ru (0001) surface. Starting from the 
optimum value of the N-Ru layer distance, we then per- 
formed an automated refinement varying the structural 
parameters mentioned above and the real part of the 
muffin-tin potential. The optimum r-factors reached for 
each of these high-symmetry configurations are compiled 
in Table I, from which it is obvious that in both phases 
the hep site represents the actual adsorption site. This 
result is consistent with a recent STM study in which 
the adsorption site of nitrogen was determined by label- 
ing hep-sites with coadsorbed oxygen atoms [22]. 

The structural parameters of the best-fit models are 
presented in Figs. 1 and 2, and the good agreement be- 
tween calculated and experimental LEED IV data is il- 
lustrated in Figs. 3 and 4. The ehemisorption of nitro- 
gen induces small lateral displacements by 0.06 ± 0.05 A 
(2 X 2) and 0.02 ± 0.05 A x \/3)R30° in the topmost 
Ru layer which are small if compared to corresponding 
displacements of about 0.1 A due to the ehemisorption 
of oxygen on Ru(OOOl) [23]. The N-Ru layer spacing 
turned out to be 1.05 A and 1.10 A in the (2 x 2) and 
the (v^ X \/3)R30° phases, respectively. These values 
are substantially smaller than those found with oxygen 
phases (about 1.24 A), which is attributed to a stronger 
bonding of N to Ru. In TDS, however, nitrogen comes 
off the surface in the temperature range from 500 K to 
850 K, while oxygen desorbs above 1000 K. Yet, both 
experimental observations are not conflicting since the 
heat of formation of molecular nitrogen (9.77 eV) is much 
higher than that of molecular oxygen (5.16 eV). Note that 
TDS yields only information about the activation energy 
of desorption which, to a first approximation, is the dif- 
ference of the Ru-X binding energy and one half of the 
binding energy of X2 (X either N or O). Ab-initio calcu- 
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lations presented here support this view and provide the 
Ru-N binding energies in the (2x2)-N, (VSx a/3)R30°-N, 
and (Ixl)-N phases to be 5.82 eV, 5.59 eV, and 4.90 eV, 
respectively, which arc indeed higher than those found in 
the (2 X 2j-0, (2xl)-0, and (Ixl)-O phases (5.55 eV, 
5.28 eV, and 4.80 eV) [24]. In all these cases, DFT cal- 
culations showed that adsorption takes place in the hep 
site. Note that the (Ixl)-N phase is not thermodynami- 
cally stable as discussed below. 

In contrast to the (2 x 2)-0 and (2xl)-0 phases on 
Ru(OOOl), which show a slight contraction of the top 
Ru-Ru layer spacing, the averaged top Ru-Ru interlayer 
spacing in the (\/3 x \/3)R3G°-N phase is slightly ex- 
panded by 1 % indicating that moderate amounts of ni- 
trogen are sufficient to deplete the population of bond- 
ing d orbitals between the first and second Ru layer, 
thus removing the contraction of the clean Ru (0001) 
surface. The hard-sphere radii of N found on the 
Ru (0001) surface are 0.58 A in both phases which com- 
pares well with values already found on other metal sur- 
faces, such as Ni(110)-(2x3)-N (0.59 A) [25], Cu(llO)- 
(2x3)-N (0.62 A) [26], and Rh (110)-(2x 1)-N (0.57 A) 
[27]. 

In Table II, the main structural parameters of both ni- 
trogen phases obtained by ab-initio calculations are com- 
pared to those obtained by LEED. Within the quoted 
error bars these parameters are identical. Even details 
of the N structures are reconciled with DFT calculations 
such as the atomic displacements of the Ru atoms in 
the (2 X 2)-N phase: While LEED revealed a rumpling 
of 0.07 ± 0.04 A and a lateral shift of 0.06 ± 0.05 A, 
DFT calculations determined these displacements to be 
0.086 A (vertical) and 0.067 A (lateral). In addition to 
these experimentally observed N surface structures, DFT 
calculations for a hypothetical (Ixl)-N structure were 
carried out. The Ru-N layer spacing (1.24 A) turned 
out to be larger than that for the (2 x 2)-N and the 
(\/3 X \/3)R30°-N phases, and the Ru-Ru layer spacing is 
heavily expanded by 5.2 % with respect to the bulk value. 
This investigation was triggered by a recent DFT /LEED 
investigation of the (Ixl)-O phase on Ru(OOOl) [24]. 
Also with the (Ixl)-O structure it was found that the 
first Ru Ru layer distance is substantially expanded by 
3.5 %. The main difference between the (Ix 1)-N and the 
(Ixl)-O phases on Ru (0001) is, however, that according 
to DFT calculations, the (Ixl)-O is exothermic with re- 
spect to molecular oxygen, while the (Ixl)-N phase is 
not. The binding energy of nitrogen in the (1x1) phase 
is only 4.52 eV which is by 0.6 eV smaller than half the 
formation energy of molecular nitrogen. 

To determine the critical N-covcragc which can be ther- 
mally stabilized on the Ru (0001) surface, we performed 
DFT calculations of the (2xl)-N and the (2x2)-3N over- 



layers assuming hep- adsorption only. It turned out that 
nitrogen in the (2xl)-N is stabilized by 0.2 eV with re- 
spect to half the formation energy of N2, while nitrogen in 
the (2x2)-3N phase is endothermic by 0.2 eV. This result 
suggests the nitrogen saturation coverage to be about 0.5 
(so far not observed in experiments) which is also com- 
patible with the formation of heavy wall domain bound- 
aries as proposed for N overlayers with coverages larger 
than 0.33 [8]. 

The strong bonding of nitrogen to Ru (0001) [exem- 
plified for the (2x2)-N overlayer] can be explained most 
easily by using a simple two-level tight-binding model 
[28]. The same model has recently shown to work for 
O/Rh(110) [29] and O/Ru(0001) [30]. In this model the 
2p orbitals of N interact with d-orbitals of Ru close to 
the Fermi energy (cf. Fig. 5) , forming a bonding hybridc 
at an energy slightly lower than of the 2p orbital of the 
free N atom and an anti-bonding hybride whose energy is 
3 eV above the Fermi energy. The bonding hybride con- 
sists (mainly) of sp hybride orbitals of the N atom and 
dz^-like orbitals of the three nearest-neighbor Ru atoms, 
aligned along the bond axis between the N and Ru atoms. 
In the non-bonding hybride, at energies — 3.5 eV < e < 
—2.5 eV below Fermi, there are one-electron states with 
pz-character at the N atom and d character at the Ru 
atoms. Close to the Fermi energy many d2^-like orbitals 
have been emptied upon forming the anti-bonding hy- 
bride, at 3 eV above Fermi. The corresponding wave 
functions exhibit a distinct node between the N atoms 
and the nearest-neighbor Ru atoms thus indeed demon- 
strating anti-bonding character. 

The strong bond of N on Ru (0001) relies on two effect. 
First, a binding hybride is formed, increasing the charge 
density at the N atom (as is consistent with the positive 
change of work function upon adsorption [31]). Second, 
the anti-bonding states are pushed above the Fermi level 
and hence do not weaken the adsorbate-substrate bond. 

In conclusion, the atomic geometry of nitrogen ad- 
sorbed on Ru (0001) has been explored by making use of 
LEED and DFT calculation. Both techniques revealed 
that nitrogen atoms adsorb in the threefold hep sites on 
Ru(OOOl), the Ru-N bond length in the (2 x 2)-N and 
the (\/3 X \/3)R30°-N phase being about 1.93 A. The 
small Ru-N bond length is indicative of a strong binding 
which is also evidenced by DFT calculations, indicating 
binding energies of 5.82 eV and 5.59 eV for the (2 x 2)- 
N and the (a/3 x y3)R30°-N phases, respectively. DFT 
calculations show that a hypothetical (Ixl)-N phase on 
Ru(OOOl) is not stable thermodynamically, while a hy- 
pothetical (2xl)-N phase would be stable. 
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TABLE I. Optimum Pendry r-factors for different structural models of the Ru (0001)-N-(2 x 2) and 
Ru (0001)-N-(V3 X V3)R30° 





(2 X 2)-N 


{VS X V3)R30°-N 


on-top 


0.53 


0.63 


bridge 


0.46 


0.52 


fee 


0.70 


0.59 


hep 


0.24 


0.29 
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TABLE II. Comparison of structural paxameters obtained by LEED and DFT calculations. The expansion of the first Ru 



layer spacing is given with respect to the bulk layer distance. 



method 


LEED 




DFT calculations 


phase 


(2 X 2)-N 




(2 X 2)-N 


{V3 X VS)-N 


Ru-N 


(1.93 ±0.05) A 


(1.93 ±0.06) A 


1.94 A 


1.94 A 


Ru-N layer spacing 


(1.05 ±0.05) A 


(1.10 ±0.06) A 


1.08 A 


1.12 A 


Ru expansion 


(-2 ± 2) % 


(+1 ± 2) % 


% 


-0.5 % 


adsorption site of N 


hep 


hep 


hep 


hep 
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FIG. 1. Hard-sphere model of the (2 x 2)-N structure on Ru(OOOl). Atomic nitrogen resides in the threefold hep hollow 
position. The Ru-N bond length obtained by LEED is 1.93 ±0.05 A while the value provided by DFT calculation is 1.94 A. 
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Ru(0001)-(V 3xN 3)R30P -N 




(0.02±0.05) 




{1.1010.06) 
(2.1610.05) 

2.14 
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FIG. 2. Hard-sphere model of the (\/3 x \/3)R30°-N structure on Ru(OOOl). Atomic nitrogen resides in the threefold hep 
hollow site. The average Ru-N bond length is found to be 1.93 ± 0.06 A and 1.94 A when using LEED and DFT calculations, 
respectively. 



10 



11 




Energy (eV) 



12 



FIG. 3. Comparison of experimental and calculated LEED I(V) curves of the Ru(0001)-(2 x 2)-N structure. The overall 
Pendry r-factor is Rp = 0.24. 
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FIG. 4. Compaxison of experimental and calculated LEED I(V) curves of the Ru (0001)-(a/3 x a/3)R30°-N structure. The 
overall Pendry r-factor is Rp = 0.29. 
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FIG. 5. The difference in density of states, n(Ru+N) - n(Ru), where n(Ru+N) is the density of states of the adsorbate 
system Ru(0001)-(2 x 2)-N and n(Ru) that of the clean Ru surface. The one-electron eigenvalues are broadened by 0.6 eV. 
The density of a representative bonding, non-bonding and anti-bonding orbital are shown in the insets. The arrows point to 
the one-particle energy regions for which the orbitals are shown. The positions of the N and Ru atoms are indicated by small 
and large circles, respectively. 
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